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A scheme of laser wakefield acceleration, when a relatively rare and long bunch of non-relativistic 
or weakly-relativistic electrons is initially in front of the laser pulse, is suggested and considered. The 
motion of test electrons is studied both in the one-dimensional case (ID wakefield) and in the three- 
dimensional laser wakefield excited in a plasma channel. It is shown that the bunch is trapped, effectively 
compressed both in longitudinal and transverse directions and accelerated to ultra-relativistic energies 
in the region of first accelerating maximum of the wakefield. The accelerated bunch has sizes much less 
than the plasma wavelength and relatively small energy spread. 

PACS numbers: 41.75.Jv, 52.35.Mw, 52.75.Di 



I. INTRODUCTION 

The rapid progress in the technology of high-intensity lasers, based on the chirpcd-pulsc amplification 
(CPA) [1], opens new opportunities for the use of lasers in many branches of science and industry. 
Relatively inexpensive tabletop terawatt lasers (so-called T 3 -lasers) become a qualitatively new tool 
in physical researches and now are available at many laboratories over the world. Last years CPA 
technique permitted the production of subpicosecond laser pulses of multiterawatt power with peak 
intensities exceeding 10 20 W/cm 2 [2]. With intensities as such we practically have to do with a new 
interaction range of laser radiation with matter, where the role of the nonlinear effects is often essential. 
In this intense laser field the matter is usually transformed to plasma and free electrons oscillate with 
relativistic quivering energy. Presently, the interactions of high-power laser radiation with plasma are 
actively investigated in connection with different applications: the excitation of strong plasma wake waves 
for focusing and acceleration of charged bunches [3] ; generation of radiation at harmonics of carrier laser 
frequency [4]; X-ray sources [5]; laser inertial fusion [6] etc. 

The laser wakefield, generated in plasma by the short (with the length « A p /2, where A p is the plasma 
wavelength) intense laser pulse provides the acceleration gradient up to tens GeV/m (laser wakefield 
acceleration, LWFA [7,8]), that is three orders of magnitude higher than that achieved in conventional 
accelerators. The main aim of experimental and theoretical works, that are presently in progress, is the 
construction of compact and relatively inexpensive accelerators of charged particles for applications in 
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physics research, medicine and hi-tech industry. However, some challenges remain on this way, one of 
the main of those is the problem of electron bunch injection. 

The wake wavelength in the LWFA is A p w 2ct£ [3] (here tl is the laser pulse duration) and makes 
up tens or hundreds micrometers for typical plasma densities n p ~ 10 16 — 10 19 cm~ 3 . To obtain high 
quality relativistic electron bunch accelerated by the wake wave, it is necessary to inject short (with the 
length L <C X p ), enough dense relativistic electron bunch in the accelerating phase of the wake wave 
with femtosecond synchronization, that is difficult technical problem (see e.g. Ref. [9]). The injection 
schemes proposed for the standard LWFA (the LIPA [9], the colliding laser pulses [10] and the LILAC 
[11] schemes) are aimed at the generation of such a short relativistic bunch. 

The diffraction broadening leads to rapid decrease of the intense laser pulse amplitude with the 
characteristic length Zr = tttq/Xl (here Zr is the Rayleigh length, ro is the focal spot size of the 
pulse, and Xl is the laser wavelength) that is typically in order of a millimeter. To prevent diffraction 
the plasma channel, with minimum density at the axis, proposed to guide the laser pulse in the LWFA 
[12]. The amplitude of the accelerating component of the wake wave, generated in the plasma channel, 
decreases as the distance from the laser pulse increases [13,14]. Besides, the change of the plasma 
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wavelength A p in the transverse direction [X p ~ n p (r), where n p is the density of electrons in the 
plasma channel], leads to undesirable wave front curving; this effect becomes stronger as the distance 
from the pulse increases. The effect of the wave front curving in the channel, in the case of a strong laser 
pulse (ao = eE m /m e cu>L > 1, where e and m e are the absolute charge and mass of the electron, E m is 
the maximum amplitude of laser field, c is the speed of light in vacuum, and u>l is the laser frequency) is 
amplified due to the nonlinear change of the wake wavelength in transverse direction [15,16]. Thus, for 
regular acceleration of a charged bunch in the wake wave, the most preferable is the region of the first 
maximum of accelerating field behind the laser pulse. 

To avoid the aforesaid difficulties in LWFA, we suggest and study in this work a scheme of trapping, 
compression and acceleration of a non-relativistic or weakly-relativistic electron bunch in the laser wake- 
field, when the bunch is initially in front of the laser pulse. The initial bunch density can be much less 
than that required for the accelerating bunch and the bunch sizes - in order or more than the plasma 
wavelength, i.e. much more than require other methods of injection [9-11]. Our investigations take into 
account both the pulse ponderomotive force and the wakefield. It is shown that the bunch is trapped, 
effectively compressed both in longitudinal and transverse directions and accelerated to ultra-relativistic 
energies in the region of first accelerating maximum of the wakefield. The accelerated bunch has sizes 
much less than the plasma wavelength and enough good quality. 

II. THE CASE OF WIDE LASER PULSE 

At first we neglect the transverse variation of the laser pulse amplitude and consider the case of 
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one-dimensional laser wakefield excited by a wide bunch in uniform plasma. This allows to study the 
longitudinal dynamics of the bunch electrons in more details. 

A. Basic equations and correlations 

One-dimensional steady wakefield excited by the linearly-polarized laser pulses are described by the 
following equation (see, e.g. Ref. [3]) 

d 2 *, R .L a ^/(l + a 2 ^) 1 / 2 ) 

where $ = 1 + ef/m e c 2 is the dimensionless potential of the plasma wakefield, a = eE {^)/m e ojjL, E 
is the electric field amplitude of the laser pulse, £ = k p (z — v g t), k p = uj p /v g , uj p = (47rn p e 2 /m e ) 1 / 2 is the 
plasma frequency, v g is the group velocity of the laser pulse which is equal to the phase velocity of the 
wake wave, fi g — v g /c, 7 = (1 — /3 2 )~ 1//2 is the relativistic factor, which, in the case 7^ >> 1, is nearly 
equal to uol/uj p . The electric field of exited wakefield, normalized to the non-relativistic wave-breaking 
field Ewb — m e v g ujp/e, can be obtained from equation E z = — (l//3 g ) 2 <i<I>/d£. The equation of motion 
of the test electron in the wakefield and in the field of laser pulse is (sec e.g. Ref. [f 7]): 

dp _ 1 da 2 

Here the first term on the right-hand side is the relativistic ponderomotive force averaged over the fast 
laser oscillations, and the second one corresponds to the plasma wakefield excited by the laser pulse, 
P = v/c, p = /?7 and 7 = (1 + p 2 + a 2 /2) 1 / 2 = [(1 + a 2 /2)/(l - /3 2 )] 1 / 2 are the normalized longitudinal 
velocity and momentum and the relativistic factor of the test electron (transverse velocity is zero in this 
Section), t — uj p t. Multiplying Eq. (2) by /?, after some simple transformations, one can obtain the 
following integral of motion (see also Refs. [14,18]) 

7 - (3 g p - $ = const. (3) 

Let us consider an electron that is initially ahead of the laser pulse at some point £0 where $ = 1 and 
a = 0. If the electron has initial momentum p < (3 g j g , it will be overtaken by the laser pulse and can 
be trapped at some point £ r inside the pulse or in the wake and accelerated. At the trapping point (or, 
in other words, at the point of reflection) the velocity of the electron become equal to v g . Then from 
expression (3) we have: 

S = (1 + a 2 /2) 1 /V 7s - ($ r - 1) = (1 + P 2 ) 1/2 - P 9 P- (4) 

In Eq. (4) a r and <Ev are the amplitude of the laser pulse and the wakefield potential at the reflection 
point £ r . From Eq. (4) it follows 
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P = l 9 [P 9 l g S±{l 2 g S 2 -l) 1 l 2 ]. (5) 

The minus sign in (5) corresponds to the initial momentum po of an electron which has momentum (3 g j g 
at the point £ r and the plus sign corresponds to the momentum of free electron which initially was at 
the point £ r . Expression (5) describes both trapped and passing particles. In the wake, electrons can be 
trapped only in the region where E z < 0. Equation of motion (2) can be rewritten in the form 

where £ is the coordinate of a test electron in the frame commoving with the laser pulse. The dimen- 
sionless velocity of the electron one can obtain from expression (i = /3 S (1 + d^/dr). 

B. Numerical results 

Eqs. (1) and (6) were solved numerically for the Gaussian laser pulse, 

a = a exp[-(£ - £ c ) 2 /<t 2 ]. 

In Fig. 1 a laser pulse with clq = 2 and nonlinear wakefield excited by it are presented (here and below 
in numerical calculations a z — 2, £ c = 3 and j g = 10). The amplitude of the wake wave is essentially 
less than one-dimensional relativistic wave-breaking field E re i = [2(1 — 7 g )] 1,/2 //3g ~ 4.26 [3,19]. Figure 2 
shows the dependence of initial electron momentum p on the trapping point near the first accelerating 
maximum in the wake wave. The minimum value of the initial momentum p m ; n corresponds to the 
trapping point where the potential is at the minimum and E z = 0. Curves 1 and 2 in Fig. 2 reach 
their minimums at different points, that is the consequence of the nonlinear increase of wake wavelength 
with the amplitude (this dependence can be found in Rcf. [20]). The curves were obtained numerically 
and coincide with the expression (5) for the trapped particles. Fig. 3 shows the dependence of the 
value of Pmin and wake wave amplitude E Ztinax on ao- One can see that a laser pulse with ao ~ 1 (that 
corresponds to the peak intensity of the pulse I ma x ~ 10 1S W/ cm 2 for Az, = l//m, and / max ~ 10 16 W/ cm 2 
when Ax = 10/ito) provides trapping of initially non-relativistic or weakly-relativistic electrons in the 
wake wave. For example, p m - m ~ 0.4 for the wakefield presented in Fig. 1. Electrons with p < p m - m 
can not be trapped in the wake wave and may be detected behind the wave. This circumstance can 
help to determine the wake wave amplitude in experiments. Our numerical calculations have witnessed 
that electrons with pa < P g "f g can not be trapped in the region occupied by the laser pulse because of 
the decelerating wakefield; only electrons with po ~ flglg are trapped in the head of the pulse (where 
E z w 0) due to the ponderomotive force. This confirms with the results of Rcf. [18]. 

Figure 4 shows the behavior of electrons of mono-energetic non-relativistic (po — 0.5, 70 = (1 + 
Pp) 1 / 2 « 1.12) bunch in the wakefield presented in Fig. 1. Initial dimensionless bunch length is L Q = 5, 
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that roughly corresponds to the linear plasma wavelength X p . When r = 50, the trapped bunch length 
is L w 0.027 and L w 0.04 when r — 100, that is two orders of magnitude less than the initial bunch 
length. The absolute energy spread A7 in the accelerating bunch increases insignificantly with time, but 
the relative energy spread e = A7/7 falls due to growing 7; for example e « 0.26 at r = 50, and e w O.f 4 
when r = 100. The acceleration gradient in the considering case is approximately equal to 2MeV/ X p . 
For example, when X p = 100/xm (n p « 10 17 cto~ 3 ), the acceleration gradient is 20GeV/m. 

Figure 5 shows the motion of electrons with different initial momcntums and the same initial positions 
(0.6 < po < f .2, 1.17 < 7 < f .56, £0 = 0) in the wakefield presented in Fig. 1. The trapped bunch length 
is nearly 27 times less than the plasma wavelength X p . The relative energy spread at r = fOO is about 
O.f, that is much less than that of initial electrons. 

The dephasing length, for electrons with p m i n < po < 1 .2, varies in the range 630 < Ld < 700 (the 
grater values correspond to the smaller initial momentums) that is comparable with the linear dephasing 
length A P 7g [3], which, in our notations, corresponds to Ld = 2^7^ = 2007T. The maximum relativistic 
factor of accelerated particles is in the range 350 < 7 max < 410 (here again the greater values correspond 
to smaller po) that essentially exceeds the linear value 27^ = 200 [3], but is an order of magnitude less 
than the maximum nonlinear value 47^ = 4000 [20,21]. 

C. Energy spread in the accelerating bunch 

Energy spread in the trapped bunch depends on energy spread and length of the initial bunch. The 
tail electrons of initial bunch are trapped earlier and therefore, have greater energy during acceleration 
(see Fig. 4). Slower particles also are trapped earlier (see Fig. 5). Let us suppose that initially the 
bunch is at the head of laser pulse, so that £ = corresponds to the bunch tail, and T tr (po) is the 
time necessary to trap an electron which is initially at £ = 0; the initial electron momentum is in the 
range pi < p < p 2 - Then, for energy spread in the trapped bunch we can write A7 ~ AT tr E z mlix = 
[Ttripz) — Ttr(pi) + Lq / (1 — V2 /v g )] £ z .max, where Artr is the time interval which is necessary to trap the 
initial bunch. For the relative energy spread one has: e ~ Arj r /(r — Art r ). These estimates agree well 
with the numerical results. One can see that the presence of fast electrons (with vo ~ v g w c) in the 
initial bunch leads to an undesirable increase in the energy spread. 

The trapped bunch density can be found from expression nfe(r) « u^Lq / 'L(t) , where n^o is the initial 
bunch density. 

D. Wakefield generated by the accelerating bunch 

The trapped bunch also generates wakefield which can destroy the laser wakefield and decrease the 
accelerating field. Because the accelerating bunch is short (L(t) <C A p ) we can consider it as a plane 
bunch and find the normalized amplitude of the wakefield excited by the bunch from expression i?b. max = 
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k p (vb/c)(Nb/n p ) [22], where Vb and Nb are velocity and the surface density of the bunch correspondingly. 
This expression is valid both in linear and non-linear regimes. In our case Nb = 6riboLo/k p , where 5 < 1 
is the ratio of number of trapped electrons to the total number of particles in the initial bunch, and we 
have: 



£&,max = S(v b /c)(n b0 L /n p ). (7) 

The normalized amplitude of moderately nonlinear laser wake wave, considering in this paper, is about 
unit. So, we can neglect the wakeficld generated by the bunch if Eb, max <C 1, or when 

n bQ < n p (c/v b ){l/5L ). 

For n p ~ 10 16 — 10 18 cto~ 3 (that is typical for the LWFA experiments [3]), Vb « c, S w 1 and the initial 
bunch length in order of X p (L ~ 5 — 10) this condition reads ribo < 10 14 — 10 16 cm~ 3 . The density of 
accelerating bunch may be in order of plasma density. 

Thus, the one-dimensional analysis has showed the possibility of trapping, essential compression and 
high-gradient acceleration of a low energy electron bunch in moderately nonlinear laser wakefield. 

III. TRAPPING, COMPRESSION AND ACCELERATION IN THE LASER 
WAKEFIELD EXCITED IN PLASMA CHANNEL 

In this section we consider our scheme of LWFA for the case of laser wakefield excited in a plasma 
channel and study the peculiarities of radial motion of test electrons during trapping and acceleration. 

A. Nonlinear laser wakefield excited in plasma channel 

As was mentioned in Introduction, the plasma channel is necessary to guide a laser pulse. This 
allows to essentially increase the laser-plasma interaction distance [12], that, in its turn, provides ultra- 
relativistic acceleration in the wakefield [3]. Nonlinear axially-symmetrical laser wakefields excited in a 
plasma channel are described by the following system of equations [15]: 

4^-^-/^ = 0, (8.2) 

'^+^+^ = 0, (8.3) 
BE 

V^H e +[3^+f3 z N e = 0, (8.4) 
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dH e dE dE z 
""flf " "flf + = °' (8 ' 5) 

N e = N p (r) - V±E r - (8.6) 

where E z ^ r and Hg are longitudinal and radial components of the electric field and azimuthal component 
of the magnetic field normalized to the on-axis wave-breaking field EwB{r = 0) = m e Lu p (r = 0)v g /e, 
p z . r are the normalized components of plasma electron momentum, 7,= = (1 + p z + p 2 + a 2 /2) 1 / 2 is the 
relativistic factor, /3 Zi r — Pz. r lie, N e = n e (£, r)/n p (0) is the normalized density of plasma electrons, 
n p (r) is unperturbed plasma density in the channel, N p — n p (r)/n p (0), Vj_ = d/dr + 1/r. The force 
acting on the relativistic electrons in the wakefield is F(— eE z , — e(E r — /3H#), 0). According to (8.5) 

dE z = d(E r - f3H e ) = _dfr 

dr d£ ~ d£ ' U 

So, the field of forces F is potential because y x F = , and one can write F = y<&(£,r), here 
<P = I - f° E z dt 

In this section we consider an axially-symmetric laser pulse which has Gaussian profile both in 
longitudinal and radial directions: 

a(£, r) = a exp[-(^ - £c) 2 /v 2 z ] exp(-r 2 /cr, 2 ). 

The laser pulse is guided in preformed plasma channel which has the following unperturbed electron 
density: 

N p= ( 1 + A ^) eXP (^)' (10) 
where r c h, A and i< 1 are constant values. Such a density profile is typical for plasma channels created 
in experiments [23]. Suppose that the pulse is guiding without change in its radius a r . In this case 
c r = T c h and n p (r c h) — n p (0) = l/irr e rl h , where r e — e 2 /m e c 2 « 2.8 x 10 _13 cm is the classical electron 
radius and all values are dimensional [12]. Then, in expression (10), A = (2/a r f3 g ) 2 . 

Equations (8.1)-(8.6) were solved numerically for the following parameters of the problem: ao = 2, 
<r z = 2, ay = 5 and -y g = 10. In this case A w 0.16, the value of b was chosen to be 0.01. In Fig. 7 we 
present the radial profile of unperturbed plasma density and the radial behavior of laser pulse intensity, 
namely exp(—2r 2 /a 2 ). Fig. 8 shows the longitudinal electric field and the focusing field f r — 0H$ — E r 
of the wakefield excited. One can see that the wake wavelength decreases as r increases. This is caused 
by the radial increase of unperturbed plasma density in the channel [3,13] and by the nonlinear increase 
of wavelength with the wake wave amplitude which is at maximum on the axis [15,16,24]. Fig. 7 shows 
also the nonlinear steepening of the accelerating field like that takes place in one-dimensional wakefield 
(see Fig. 1). Due to the dependence of the wavelength on r, the field in the radial direction grows more 
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chaotic as the distance from the laser pulse increases. In fact, the oscillations of the plasma for different 
r are started behind the pulse with nearly equal phases but different wavelengths. As |£| increases, the 
change of phase in the transverse direction becomes more and more marked. This leads to a curving of 
the phase front and to oscillations in the transverse direction [15,16,24]. Such behavior of the wakefield 
excited in a plasma channel leads to the transverse multistream motion of plasma electrons in the wake 
and to the transverse wave-breaking [25]. The radial dependence of longitudinal electric field and the 
focusing force is shown in Fig. 8 for point £ = —10.9 at which the on-axis accelerating field reaches 
its maximum. We see that the wakefield changes its sign and is steepened. For the ultra-relativistic 
acceleration of electrons one needs to use a region in the wakefield where the conditions E z < and 
f r < are satisfied simultaneously. The radial steepening leads to the radial restriction or the region 
suitable for acceleration. Near the first accelerating maximum of the wakefield shown in Fig. 7, the 
suitable region is r < 2.8. As the distance from the laser pulse increases, the suitable region becomes 
narrower, so that at some distance the wakefield is highly irregular. Thus, the most preferable for electron 
acceleration is the region of the first accelerating maximum in the wake. 

B. Equation of motion of bunch electrons 

Three-dimensional vector equation of motion of bunch electrons is 

Here all values are dimensionless, (3 = v/c = p/7 is the normalized velocity, 7 = (1 + p 2 + a 2 /2) 1 / 2 = 
[(1 + a 2 /2)/(l — /3 2 )] 1 / 2 is the relativistic factor. For the momentum components, from Eq. (11) one 
has: 



*»=0, (12.2) 
dr 



t=-u E > + ^~wM' (i2 - 3) 

It follows from Eq. (12.2) that the azimuthal momentum is conserved, pe — const, /3e(r) = ps(0)/j(t). 
The azimuthal momentum has not essential influence on the axial and radial dynamics, and we suppose 
Pe(0) — in this paper. For the energy of electrons Eq. (11) gives the following equation: 
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From Eqs. (12.4), (13) and (9) we obtain the integral of motion 



7 - f3 gPz - r) = const, (14) 

which formally coincides with the one-dimensional integral of motion (3) [14,26,27]. Electrons can be 
trapped in the region where wakefield is both accelerating and focusing. For the scattered particles, from 
Eq. (14) one has p 2 = (S + /3 g p z ) 2 - p 2 z - 1, here S = [1 + p 2 (0)] 1 / 2 - ( 3 g p z (0). If an electron is initially 
non-relativistic (|p(0)| < 1, S « 1), then p r w (2p z ) 1 / 2 and tan6> = p r /p z « [2/(7 - 1)] 1/2 , where 6> is 
the angle between z-axis and final momentum of the scattered electron [28] . 
Taking into account Eq. (13), we rewrite Eqs. (12.1) and (12.3) in the form 

+ i [(1 - 0>)E, - ,W6 + 3 r H,\ + 1L =g|^ 1 ff = 0. (15.1) 

% + I [(1 - * )* - AAA - &«] - ^ ( A g - . 0. ,15.2) 

The normalized components of velocity obey the formulae (3 Z = + d^/dr) and /3 r = P g dr/dr. For 
particles trapped in the wakefield, we suppose that during acceleration /3 Z w 1, /3 2 <iC 1 and r < 1 (the 
numerical results presented below show that this is the case). Then, from Eq. (15.1) one has 

d 2 ^/dT 2 « i? z /7 3 . (16) 

It follows from this equation that d'y/dr w — i? z and 7 ~ ~ f E z dr. Thus, the longitudinal dynamics of 
accelerating particles is approximately the same as in the one-dimensional case. The radial motion of 
electrons, according to Eq. (15.2), obeys the equation 

£ + S* +f! ^0, (17, 
dr z 7 dr 

where Q. = (|9/ r /9r|/7) 1 / 2 is the betatron frequency. Supposing that the value of E z is approximately 
conserved during acceleration, we can write 7 w |^|(t — r tr ). In this case solution of Eq. (17) is 

r = r(T tr )J [2(\df r /dr\(r - T tr )/\E Z \)V% (18) 
where Jo is the Bessel function of zero order. 

C. Results of test-particle simulations and discussion 

Motion of test electrons in the 2D wakefield presented in Fig. 7 was investigated by numerical solution 
of Eqs. (15.1) and (15.2) for different initial positions and momcntums. Figure 9 shows the behavior 
of electrons with zero initial transverse momcntums and with different initial radial positions. One can 
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see that particles are trapped near the first accelerating maximum in the wake. During the trapping, 
electrons concentrate near the axis due to the focusing force (3 z Hg — E r . Sins the longitudinal size of 
the trapped bunch is much less than the plasma wavelength and its transverse size is essentially less 
than that of the laser pulse, the electrons experience approximately the same accelerating field; the 
longitudinal dynamics is well described by the one-dimensional theory. The focusing force acting on 
the bunch electrons depends on r linearly (see Fig. 8). The small bunch sizes (as compared with the 
wakcfield characteristic sizes) and the fact that electrons are trapped near the accelerating maximum 
provide high accelerating gradient and relatively small energy spread. For example, the relative energy 
spread of electrons presented in Fig. 9 is 5% at r = 300. The numerical results show that dynamics of 
the accelerating bunch is well described by approximate equations (16)-(18). The betatron oscillations 
of the accelerating electrons are clearly seen in Fig. 9(b). The wavelength of this oscillations decreases 
with the increase of particle's energy that conform to the formula for betatron frequency. Radial velocity 
of accelerating electrons is much less than the longitudinal one, |/3 r (r)| < 0.1. One can see also that 
even electrons which are initially at the periphery [r(r = 0) = ro ~ ay] can be trapped in the wakefield 
and accelerated. The characteristic dependence of the niininium trapping threshold Pz,m'm 

on the initial 

radial position of an electron is presented in Fig. 10. Figure 11 shows the minimum and maximum initial 
radial momcntums of trapped electrons in dependence on initial radial position. The figure witnesses 
that electrons which initially move at relatively high angle to the axis (up to tens degrees) also can 
be trapped and accelerated. This again is caused by the focusing force and the fact that electrons are 
initially non-relativistic (70 ~ 1). 

In Fig. 12 we show behavior of electrons of a bunch with the following initial parameters: < £0 < 5, 
r o < 4, 0.6 < p z o < 0.8, —0.02 < p r o < 0.02. The passing particles (not showed) are well separated from 
accelerating one both spatially and energetically. The length of accelerating bunch in this case also is 
much less than the plasma wavelength [L(t = 100) ss 0.27, L(t = 300) « 0.19]. The radius of the bunch 
R decreases relatively slowly during acceleration and is essentially less than the characteristic transverse 
size of the wakefield ay, R(t) ~ 1; the bunch radius can be reduced by the choice of smaller laser spot 
size. The absolute energy spread does not change practically, A7 « 24, but the relative energy spread 
falls and is equal to about 10% at r = 300. The estimations of absolute and relative energy spreads 
presented in Sec. II are valid also in 3D case. 

Total number of electrons trapped and their density can be estimated from expressions 

N tot ~ 5im W (j 2 r L Q /kl, (19) 
n b ~6n w (a r /R) 2 (L Q /L). (20) 
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The on-axis amplitude of the linear wake wave excited by the bunch is reduced by the factor T(i?) = 
1 — RKi (R) < 1 [29] (where K\ is the modified Bessel function) as compared to the one-dimensional case 
(see Sec. II). Therefore, in our case, for the amplitude of wakefield generated by the bunch, we have: 
Eb,Tnax ~ TLnb/rip. This wakefield can be neglected when -E& imax <C £ z . ma x, or taking into account (20) 
- if STL (a r /R) 2 (n b0 /n p ) <C 1; when R <C 1, T w R 2 /2, and this condition reads SL a 2 nb /2n p <C 1. 
Total number of bunch electrons, according to (19), is restricted by the following condition: N tot <C 
irnpkp 3 (R 2 /T) w 1.4 x 10 7 (R 2 /T)X p \pm]. 

For the normalized emittance e n — o 2 j(5 (here cto is the matched transverse size of the bunch, (3 is 
the betatron length) of the accelerating bunch, in our notations, one can write e n ~ R 2 Q,X p /A-k 2 . In the 
case X p = 100/iTO (n p w 10 17 cto~ 3 ), for the bunch presented in Fig. 12, e n ~ 877,771/7 1/7 2 ; for example, 
e n ~ 0.5nm when 7 = 300, that is comparable with the emittance expected in the TeV-range laser 
wakefield accelerator [30,31] (see also Refs. [14,31-34] for the dynamics of accelerating bunch). 

IV. SUMMARY 

The results of the present work show the possibility of trapping, essential compression both in lon- 
gitudinal and transverse directions and ultra-relativistic acceleration of an initially non-relativistic or 
weakly-relativistic electron bunch in moderately nonlinear (a ~ 1, E z ^ max ~ 1) laser wakefield. The 
initial bunch can be generated, for example, by a photocathode. So far as electron bunch is initially 
non-relativistic (70 ~ 1), trapping and compression take place during time interval comparable with the 
plasma wave period, that is much less than the time scale of longitudinal dynamics of relativistic particles 
in the wake [33] . Due to the fact that trapped bunch sizes are essentially less than characteristic spatial 
scales of the wake wave, the energy spread in the accelerated bunch can be relatively low, namely few 
percent. In our scheme the problems connected with the wake wavefront curvature also are removed. 
The electron bunch trapped and accelerated can be accelerated further in the multi-stage LWFA [31]. 

Thus, the scheme of LWFA proposed, has the following advantages: (a) instead of injection of an 
enough dense relativistic electron bunch with small sizes (in order of a micrometer), our scheme utilizes 
a non-relativistic, rare and long electron bunch, that is much easier to get technically, (b) femtosecond 
electron bunch synchronization in the laser wakefield is not required, (c) effective electron bunch com- 
pression, and (d) spatial and energetic separation of the initial electron bunch, that can decrease the 
trapped bunch emittance. 
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FIGURE CAPTIONS 



Fig. 1. The one-dimensional nonlinear wakcficld excited by the linearly-polarized laser pulse with 
peak normalized amplitude ao = 2, <j z = 2, 7 = 10. 1 - The normalized longitudinal electric field 

2 - the dimensionless potential of the wakefield <&(£); 3 - the amplitude of the laser pulse a(£). 

Fig. 2. Dependence of electron initial momentum pq on the trapping point near the first accelerating 
maximum. 1 - ao = 2; 2 - ao = 3. 

Fig. 3. The minimum momentum of the trapped electrons p m i n (curve 1) and the wake wave 
amplitude -E Zimax (curve 2) in dependence on peak amplitude of the laser pulse ao. 

Fig. 4. Trapping, compression and acceleration of initially mono-energetic electron bunch in the 
wakcficld presented in Fig. 1, po = 0.5, 1 < £0 < 6. Evolution of the coordinate (a) and the relativistic 
factor (b) of electrons. 

Fig. 5. Behavior of electrons with initial position £0 — and with initial momentums po = 0.5, 0.8, 1 
and 1.2 in the wakcficld shown in Fig. 1. Electrons with smaller initial momentums are trapped earlier, 
(a) Coordinate and (b) relativistic factor of electrons. 

Fig. 6. The radial profiles of unperturbed electron density in the plasma channel (curve 1) and the 
laser pulse (curve 2), r c h = cr r = 5, b = 0.01. 

Fig. 7. The two-dimensional nonlinear laser wakefield excited in the plasma channel with the radial 
density profile shown in Fig. 6, ao = 2, a z = 2, ay = 5. (a) The longitudinal electric field for tq = 0, 3 
and 5 in the order of magnitude reduction, (b) The focusing field f r = (i g Hg — E r . l-r = l;2-r = 3; 

3 - r = 5. 

Fig. 8. The radial behavior of the wakefield shown in Fig. 7, at £ = —10.9. 1 - longitudinal electric 
field E z (£, = -10.9, r); 2 - the focusing force f r (£ = -10.9, r). 

Fig. 9. Trapping and acceleration of electrons with zero initial momentums in the wakefield given in 
Fig. 7, p z o = 0.8, £0 = 0. Longitudinal (a) and radial (b) positions and relativistic factor (c) of electrons. 

Fig. 10. The characteristic dependence of the minimum trapping threshold on initial radial position 
of electron, p r0 = 0, £0 = 0. 

Fig. 11. Maximum (curve 1) and minimum (curve 2) initial radial momentums of trapped electrons 
depending on initial radial position, p Z Q = 0.8, £0 = 0. 

Fig. 12. Trapping, compression and acceleration of an electron bunch in the wakefield presented in 
Fig. 7. Initial parameters of the bunch are: < £0 < 5, ro < 4, 0.6 < p z o < 0.8, —0.02 < p r o < 0.02. 
Radial positions (a) and relativistic factor (b) of electrons. 
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